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Abstract. - We have studied the diffusion of a single particle on a one-dimensional lattice. It is 
shown that, for a self-similar distribution of hopping rates, the time dependence of the mean-square 
displacement follows an anomalous power law modulated by logarithmic periodic oscillations. The 
origin of this modulation is traced to the dependence on the length of the diffusion coefficient. 
Both the random walk exponent and the period of the modulation are analytically calculated and 
confirmed by Monte Carlo simulations. 
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Brownian Motion is a well-known phenomenon, and 
since its theoretical foundations were laid, more than one 
hundred years ago [1,2], diffusion processes and random 
walk models have been attracting the attention of re- 
searchers. The importance of random walk (RW) resides 
in the fact that it is the simplest realisation of Brownian 
Motion, with applications in almost every field of science 
where stochastic dynamics play a role [3]. It is worth to 
I'emark that, even though a random walker evolves accord- 
ing to simple rules, a considerable effort may be needed to 
Solve the dynamic problem in detail and unexpected be- 
haviours may emerge. Thus, for example, our understand- 
ihg of the mechanisms responsible for anomalous diffusion 
has been strongly fluenced by the large amount of work de- 
voted to the study of RW in non-Euclidean media, during 
the last three decades [4-6]. 

In the last years, it has been often reported that, 
sometimes, the time behaviour of a RW is modulated 
by logarithmic-periodic oscillations. These fluctuations 
has been rigorously studied by mathematicians on spe- 
cial kinds of graphs. A proof of the fluctuating behaviour 
of the n-step probabilities for a simple RW on a Sierpihski 
graph was given in ref . [7] and a generalisation to the broad 
class of symmetrically self-similar graphs can be found 
in ref. [8]. Within the physical community, it has been 
shown that, on Sierpihski gaskets, the mean number of 
distinct sites visited at time thy N noninteracting random 
walkers presents an oscillatory behaviour [9] and, more re- 



cently, detailed studies of the log-periodic modulations on 
fractals with finite ramification order, were presented in 
refs. [10,11]. 

Log-periodic modulations are not restricted to random 
walks. It is in general believed that they appear because 
of an inherent self-similarity [12], responsible for a discrete 
scale invariance (DSI) [13]. Examples of these oscillations 
have been detected in earthquakes [14,15], escape prob- 
abilities in chaotic maps close to crisis [16], biased dif- 
fusion [17,18], kinetic and dynamic processes on random 
quenched and fractal media [19-22], and stock markets 
near a financial crash [23-26] . 

In this work we analyse a minimal model of RW, which 
results in log-periodic modulations of some observables. 
The main objective is to investigate the underlying physics 
of the oscillatory behaviour mentioned above. Sometimes, 
physical phenomena can be more easily grasped with the 
help of simple models. Thus, the present study may be 
useful to determine the mechanisms involved in this kind 
of oscillations. 

For the sake of simplicity, we consider the problem of a 
single particle moving on a one-dimensional lattice. At ev- 
ery time step, the particle can hop to a nearest-neighbour 
(NN) lattice site with a probability per unit time, which 
depends on the involved (initial and final) sites only. An 
additional condition is that the forward and backward 
hopping rates between a given pair of NN sites must be 
identical. One of the basic ingredients of the model is a 
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self-similar distribution of hopping rates. As we will show, 
this leads to a length-scale dependence of the diffusion co- 
efficient that serve to explain both, the anomalous overall 
behaviour and the modulations observed in higher resolu- 
tion measurements. 

Before going into the details of the self-similar model, 
let us consider a one-dimensional periodic lattice with M 
sites per unit cell. The hopping rules stated above can 
be schematically represented by a set of barriers of height 
= c/k^'^ (ieZ), where c is an arbitrary constant and 
the probability of a hop from site j to site j + 1 (and 
from site j + 1 to site j) per unit time. An sketch of this 
structure, with A/ ~ 5, is shown in fig. 1. 
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Fig. 1: A one-dimensional periodic lattice with AI — 5 sites per 
unit cell. A schematic barrier of height c/k^'^^ is drawn between 
every pair of NN sites (i and i + 1), where fc''' is the hopping 
rate from site i to site i -I- 1 (and from site i -I- 1 to site i). The 
periodicity conditions = fc^'' are satisfied. 

The effective diffusion constant Doff for a periodic model 
can be computed following the steady-state method intro- 
duced in rcf. [27]. If a is the distance between NN sites, it 
is found that 
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(in what follows we use a = 1). The meaning of eq. (1) is 
that, for times t ^ tesc, where tesc is the average time for 
the particle to escape from the initial unit cell, the parti- 
cle mean-square displacement (A^x) satisfies the normal 
diffusion equation 



{A^x{t)) = {[x{t) ~ x(0)]2) = 2D,st 



(2) 



This is a very simple example of a length-scale dependence. 
There is an asymptotic regime, described by eq. (2), for 
\/ (A^x) greater than M, and a transitory one, not de- 
scribed here, for ^y{A^x) smaller than M. 

In order to obtain the oscillations we are interested in, 
we proceed now to modify the periodic lattice to allow for 



a greater number of kinetic regimes. Let us first introduce 
a parameter L, which is an odd natural number greater 
than 1. The model is built in stages and the result of 
every stage is called a generation. The building process is 
illustrated in fig. 2, for L = 5. 
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Fig. 2; Construction of a lattice with a self-similar distribution 
of hopping rates (L = 5) . Top: The zeroth generation (n = 0) 
with all the hopping rates equal to go. Center: The first gener- 
ation {n — 1), where the hopping rates qi appear. The corre- 
sponding barriers are separated by a distance L. Bottom: The 
second generation (n = 2) shows the emergence of the hopping 
rates 52 separated by a distance L^. In the limit n ^ 00 a 
self-similar distribution of hopping rates is obtained. For more 
details, see the text. 

The zeroth-generation lattice corresponds to the situa- 
tion in which all the hopping rates are identical (fc'*' = qq, 
yi £ Z). In this stage, eq. (2), is vahd for a/ (A^x) greater 
than 1 with Dgff = D^^^ (= go, see eq. (1)). 

In the first generation, the hopping rate k^^^ is set to 
Qi (< Qo) for every j = pL — {L + l)/2, with p integer. 
All the other hopping rates remain as in the generation 
zero. Eq. (2) is also valid in first-generation model but for 
\/ (A^x) greater than L and a new value of DcB = -D'^^ 
(also given by eq. (1)). 

This process is iterated indefinitely and, in general, the 
generation n is obtained from the generation n ~ 1 af- 
ter replacing by g„ (< In-i) the values of k'^^\ for every 
j = pi" - (L" + l)/2, with p integer. 

In the limit of an infinite number of iterations we get the 
model, with a self-similar distribution of hopping rates and 
an infinite set of diffusion constants {D^^\n = 0, 1, 2, ...}, 
discussed in the rest of the work. 
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It is not hard to convince oneself that, in th is mod el, 
eq. (2), with Deff = D^"^ , should hold for ^(A^) in 
some range between L" and and this for every non- 

negative integer n. If the RW begins from the site j = 
(the centre of symmetry of the lattice), for short enough 

times (though longer than t^Pgc ~ l/'Zo) the particle be- 
haves as being in the zeroth-generation lattice. It only 
feels the action of the lowest barriers and normal diffu- 
sion, with a constant D^'^\ should be observed. However, 
when ^y (A'^x) is of the order of L (t ~ tiW), the particle 
starts to interact with the barriers of height c/qi. For even 
longer length scales, though shorter than {t < fisc), ev- 
erything happens as in the first-generation lattice, and one 
should then observe normal diffusion with a constant D^^\ 
Because of the self-similar properties of the lattice, this se- 
quence of changes continues indefinitely, and the effective 
constant I?^"', which corresponds to the normal diffusion 
in the n*^-generation lattice, should appear at a scale L" . 
If, in addition, we impose that 

, ^ 1 + A, forn = 0,l,2,... , (3) 

^(n+l) ' 5 1 5 J V / 

where A > is another parameter of the model, both 
the diffusion coefficients and the hopping rates are deter- 
mined, up to a multiplicative constant, (see eqs. (1) and 
(3)) through 

^(«)^ 90 for n = 0,1, 2, 3,... , (4) 
and the iterative relation 
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Fig. 3: (Color online) Schematic of the mean-square displace- 
ment as a function of the time, shown by the thick red curve. 
The length of the segment be is log(2D(^)) - log(2D(2)) = 
log(l + A), because of eq. (3). From the slopes {— 1) of 
the full straight lines (representing the normal diffusion be- 
haviours, (A^x) = 2D^"h), one gets that the segments ad 
and cd have the same length or, equivalently, that log r = 
logL'^ + log(l + The dashed straight line represents the 
global power law (A^a::) ^ f^", with 2^ = log L^/ log r. More 
details in the text. 



— ^ — + {1 + Xy-^XL\ for 2 = 1,2,3... . (5) 

From the discussion in the paragraph immediately pre- 
ceding eq. (3), we can anticipate that the mean-square 
displacement behaves qualitatively as in fig. 3. This is a 
sketch of (A^x) (thick curve), which, as a function oft, has 
a power-law form modulated by a log-periodic amplitude. 
That is, 

(A2x)(i)=C<2V(0, fori>4°i, (6) 

where v is the RW exponent, and f{t) a log-periodic func- 
tion, which satisfies /(tr) = f{t), with the logarithmic 
period log(r). The value of the constant C is obtained 
by asking that the log-time average of log(/) over one 
period be zero (sec, for more details, fig. 5). We also 
observe, in fig. 3, two groups of inclined straight lines. 
On the one hand, the dashed line, which corresponds 
to the (A^x) global power-law trend, and has a slope 
of 2i/. On the other hand, the solid lines, which repre- 
sent normal diffusion in each of the different generation 
lattices, and have slopes of 1. From these slopes, it is 
clear that 2r/ = log(L^)/log(T) (dashed line), and that 
log(r) = log(L^) + log(l + A) (solid line), which gives both 
T and v expressed in terms of the parameters L and A, 



T = Li/"^ = (1 + A)l2 , (7) 



logL 

Note that, since A > 0, an anomalous diffusion appears 
{v < 1/2, see eq. (8)), and that, from the sketch in fig. 3, 
we can predict that the amplitude of the modulation in- 
creases with the increase of A or i. 

Let us remark that the self-similarity in the mean-square 
displacement, schematically shown in fig. 3 and mathe- 
matically described by eq. (6), is a direct consequence of 
the set of constraints (3). Because of these relations, the 
distance between any pair of nearest solid straight lines 
is a constant and any pair of nearest equivalent points in 
the graph (like a and b) are related by the transformation 
(t^Tt , (A^x) -^T^^'iA^x)). 

Even though we have focused on the properties of the 
mean-square displacement, a similar analysis applies to 
the average number S{t) of distinct sites visited by the 
particle, after a time t. As we are working with a one- 
dimensional lattice, S ^ y^A^a?}, and it is thus expected 
that 

S^C'f'git) , (9) 
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where the exponent i] is equal to f, g{t) is a log-periodic 
function, g{Tt) = g{t) (with r given by eq. (7)) and the 
constant C" is obtained by asking that the log-time average 
of log{g) over one period be zero. 

To check the validity of the analytical predictions stated 
above, we have performed Monte Carlo (MC) simulations, 
with qo = 1/2 and a time step St = 1. Every simulation 
begins with the particle at the center of a sixth-generation 
lattice and stops after a given number of MC steps, always 
chosen to avoid that the particle reaches the highest (6*'*- 
order) barriers. 

The numerical results of the mean-square displacement 
as a function of the time is plotted in fig. 4 for two sets 
of parameters, L = 5, A = 0.2 and i = 5, A = 5. It is 
apparent in this figure that (A^a;)(t) satisfies a modulated 
power law. 
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Fig. 4: The mean-square displacement as a function of the 
time, for L = 5, X = 0.2 (squares) and L = 5, A = 5 (circles). 
The straight lines yi and 3/2 have slopes 21/, with the values of 

V (y = 0.4732 and v = 0.3212, respectively), given by eq. (8). 
These lines are drawn to guide the eyes. 

The modulations can however be better observed in 
fig. 5, where we have plotted log((A^a;)/Ct^'^) versus 
log(t), using the same data as in fig. 4. The values of 

V were computed from eq. (8) and the constants C are ap- 
propriately chosen to have the oscillations centered around 
zero. We have drawn three straight lines. The central one 
represents the equation (A^x) = Ct^^ . The other two, 
which have slopes of 2i^±0.02 and were shifted for the sake 
of clarity, can be used for error estimation. The curvilinear 
lines are of the form Asin(27rlog(i)/ log(T) -I- a), i. e., the 
first-harmonic approximation of a periodic function with 
period log(T), where A and a are fitted parameters. From 
fig. 5 it is indeed clear that the predictions in eqs. (7) and 
(8) are consistent with the numerical findings. It is impor- 
tant to emphasise that the theoretical-simulation agree- 



ment is also excellent for other set of parameters. This 
can, for instance, be seen in fig. 6, where we have plotted 
the results of simulations for L = 3, A = 2 and L = 9, 
A = 2. Regarding the way the amplitude of the oscillation 
depends on the parameters of the model, we can see in 
fig. 5 and in the inset of fig. 6 how the former increases 
with the increase of any of the latter, which confirm our 
predictions. We would like to mention that, although not 
shown here, we have also verified that the agreement be- 
tween numerical and theoretical results is as good for the 
mean number of distinct visited sites as it is for the mean- 
square displacement. 
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Fig. 5: The logarithm of the scaled mean-square displacement 
as a function of the logarithm of the time, for 
L = 5, A = 0.2 (squares) and L = 5, A = 5 (circles) The values 
of V were obtained from eq. (8) and C are appropriately cho- 
sen constants. The curvilinear lines represent first-harmonic 
approximations of the data, ylsin((27rlogf)/(logr) -I- a). The 
period r is given by eq. (7). A and a are fitted constants. 



So far, we have considered RW's that always begin from 
the centre of symmetry (central site in fig. 2). To exam- 
ine what happens when the RW begins from another site, 
the evolution of the mean-square displacement on a lat- 
tice with L = a and A = 5 is plotted, in fig. 7, for two 
different initial conditions. In one sample (squares), the 
particle was initially located at j = (the center of sym- 
metry), in the other (stars), at j = 12 (the rightmost site 
in fig. 2-bottom). As can be expected, two clearly different 
behaviours are present at short times. However, at long 
times, a superposition of the data is observed. This can 
be understood as follows. For times much longer than the 
time for the particle to escape the second-generation unit 
cell (i.e., when \/ (A^x) ^ L^) the diffusion is governed by 
the coefficients Z)^"\ with n > 2, and becomes indepen- 
dent of the initial position. Similar results were obtained 
for several initial positions and different values of L and 
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Fig. 6: (A^a:) against t for L — 3, \ = 2 (squares) and L — 9, 
\ — 2 (circles). The straight lines yi and 1/2 have slopes 21/, 
with 1/ (i^ = 0.40 and v = 0.33, respectively) obtained from 
eq. (8). The inset shows a plot of log{A^ x) /Ct^'^ vs. logt for 
the same data. The curvilinear lines were obtained in a similar 
way as in fig. 5. 

A. In general, it can be said that for ■\/ (A^x) ^ i", the 
behaviour of the RW does not dependent of the initial po- 
sition, provided that the latter is at a distance less than 
{L" — l)/2 from the centre of symmetry. 
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Fig. 7: Dependence of the mean-square displacement on the 
initial position. The squares represent the time behaviour of 
RW's that begin at site i — 0. The stars, those that begin at 
site i = 12. The inset shows log{{A^x) /Ct'^'^) versus logt for 
the same data. 

To summarise, we have presented a one-dimensional 



self-similar model of RW, which yields a log-periodic mod- 
ulated power law for the most commonly measured quan- 
tities {{A'^x){t) and S{t)). The goal of this paper it to 
introduce a minimal model that provides a simple way to 
gain insight into the effect of self-similarity on this kind of 
oscillations. As it was mentioned above, log-periodic mod- 
ulation is characteristic of fractals with finite ramification 
order [7-11]. In these objects, at a given length-scale, 
a random walker cannot easily pass from one array to 
the next, due to the small number of connections between 
equivalent structures; a difficulty that increases with the 
structure linear size (because of the larger ratio of inter- 
nal to connecting sites). In the model here, a similar effect 
occurs at a length L" , as a consequence of the decrease of 
the hopping rates at the edges of the n*'*— generation unit 
cell, when n changes to n + 1. This may be considered as 
a rough qualitative explanation of the log-periodic modu- 
lations in fractals with finite connexion order. Even when, 
for these objects, we have not calculated how the effective 
diffusion constant depends on a characteristic length, we 
expect that the knowledge of this dependence would allow 
us to obtain the values of and r, in the same way as 
for the one-dimensional model. Finally, let us stress that, 
by tuning the values of the two parameters L and A, it is 
possible to design a structure that leads to an oscillatory 
power-law behaviour with predetermined values of v and 
T (whenever ly < 1/2 and r = L^^^). So far, the parameter 
L has been restricted to odd values, because in these cases 
the lattice has a centre of symmetry. However, as fig. 7 
shows, the initial position of the random walker is not rel- 
evant for its long-time behaviour and thus the model can 
be easily generalised to include even values of L greater 
than 2. 
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